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PREFACE 


This  report  describes  work  performed  by  the  Advanced  Structural  Methods  department 
of  the  Northrop  Grumman  Corporation,  Military  Aircraft  Systems  Division,  for  the  Air 
Vehicles  Directorate  of  Air  Force  Research  Laboratory,  Air  Force  Materiel  Command, 
Wright-Patterson  Air  Force  Base,  Ohio,  under  contract  F33615-95-D-3215,  Project  No. 
0006,  “Smart  Actuation  Systems  for  Enhanced  Aircraft  Maneuver  Performance.” 

The  optimal  control  theory  is  applied  to  determine  the  solid  state  actuator  placement 
and  to  calculate  the  forces  required  to  perform  flight  maneuvers  of  flexible  aircraft  at 
desired  conditions.  The  algorithm  uses  the  Air  Force  sponsored  optimization  software, 
ASTROS,  to  generate  structural  and  aerodynamic  data.  A  few  flight  maneuver  cases 
have  been  exercised  to  show  the  merits  of  solid  state  actuators  in  6.0g  pull-up  of  a 
typical  fighter  aircraft. 

Dr.  Narendra  S.  Khot,  at  Air  Vehicles  Directorate,  Air  Force  Research  Laboratory, 
initiated  the  solid  state  actuator  feasibility  study  program  and  served  as  the  technical 
monitor.  Dr.  Kari  Appa  was  the  principal  investigator.  Mr.  John  Ausman  was 
responsible  for  implementing  the  optimal  control  algorithm  in  conjunction  with  the 
ASTROS  module  to  generate  required  data. 
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1.  INTRODUCTION 


1.1  Summary 

The  performance  characteristics  of  aircraft  largely  depend  on  the  quality  and 
distribution  of  airflow  on  the  lifting  surfaces.  By  nature,  birds  are  able  to  configure  their 
wings  in  such  a  manner  that  the  airflow  quality  and  flying  efficiency  are  at  optimum 
conditions.  Duplication  of  such  airflow  characteristics  on  manufactured  flying  vehicles 
has  been  a  desired  objective  of  airplane  designers,  beginning  with  the  Wright  brothers. 
To  simulate  bird-like  flying  characteristics,  the  lifting  surfaces  must  be  able  to  deform 
smoothly  at  appropriate  locations.  In  the  mid  1980s,  the  Air  Force  sponsored  a  mission 
adaptive  wing  (MAW)  project  to  study  aerodynamic  and  maneuver  performance 
characteristics  of  tactical  aircraft  (Refs.  1 , 2  and  3).  An  F-1 1 1  aircraft  was  selected  and 
fitted  with  hydraulic  actuators  to  deform  the  wing.  This  aircraft  was  test  flown  in  several 
mission  performances.  The  test  results  showed  overwhelming  aerodynamic 
performance  benefits  and  agility  characteristics.  However,  the  actuation  system  was 

heavy  and  expensive  to  operate  so  practical  implementation  of  this  concept  could  not 
be  realized  at  that  time. 

Recent  analytical  and  wind  tunnel  studies  sponsored  by  ARPA  and  the  Air  Force  (Ref 
4)  show  how  smoothly  contoured  control  surfaces  promote  incremental  growth  in 
suction  pressure  near  the  leading  edge.  This  has  a  beneficial  effect  on  control  surface 
effectiveness,  leading  to  enhanced  aircraft  maneuver  performance.  A  few  test  cases 
taken  from  this  study  are  presented  in  the  next  section  to  emphasize  the  need  for  active 
control  of  lifting  surface  camber  to  enhance  aircraft  performance.  This  report,  in 
subsequent  sections,  describes  an  analytical  approach,  which  can  be  used  to 
determine  optimum  wing  camber,  and  also  to  command  the  solid  state  actuators  to 
deform  the  desired  lifting  surface  at  reduced  power  requirement. 


1.2  Technical  Background 

In  conventional  aircraft,  the  leading  and  trailing  edge  control  surfaces  are  used  as 
aerodynamic  effectors  to  generate  desired  lift  distributions  on  wings  and  control 
surfaces.  The  leading  edges  are  generally  used  to  minimize  flow  separation  at 
moderately  high  angle  of  attack  flight  maneuver  cases,  while  the  trailing  edges  are 
used  to  obtain  desired  pitching  and  rolling  moments. 

Two  factors  limit  the  effectiveness  of  the  trailing  edge  control  surfaces.  First,  the  control 
effectiveness  decreases  with  increasing  dynamic  pressure  due  to  aeroelastic  effects 
(i.e.  adverse  twisting  of  the  wing).  The  second  problem  is  that  there  exists  a  massive 
flow  separation  along  the  hinge  line,  resulting  in  reduced  aerodynamic  loading  on  the 
wing  as  well  as  on  the  control  surface.  Figure  1-1  a  shows  flow  separation  due  to  rapid 
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change  in  control  surface  slope  at  the  hinge  line.  This  problem  can  be  alleviated  by  the 
use  of  hingeless  contoured  trailing  edge  surface  as  shown  in  Figure  1-1b. 


POINT  OF  FLOW 
SEPARATION 


FLOW  SEPARATION  POINT  MOVED 
FURTHER  DOWN  USING  SMART 
ACTUATOR  TECHNOLOGY 


(c)  SMART  CONTOURED  HINGELESS  FLAP  WITH  FLOW 
SEPARATION  SUPRESSION  DEVICES 


F9&CM/01 


Figure  1-1 .  Flow  Characteristics  over  Hinged  and  Hingeless  Contoured  Control 

Surfaces 


Further,  a  flow  separation  suppression  device  can  be  used  to  move  the  point  of 
separation  towards  the  trailing  edge  as  depicted  in  Figure  1-1c.  To  demonstrate  the 
difference  in  pressure  distributions  between  hinged  flap  and  smoothly  contoured  flaps, 
a  few  studies  were  conducted  using  computational  fluid  dynamics  (CFD)  method.  The 
computational  (CFD)  solutions  are  shown  in  Figure  1-2  in  which  the  solid  curve 
represents  the  pressure  distribution  computed  for  contoured  trailing  edge  surface,  while 
the  dotted  curve  denotes  the  data  obtained  using  a  hinged  flap 

Subsequently,  wind  tunnel  tests  were  conducted  in  NASA  Langley’s  Transonic  Dynamic 
Tunnel  (TDT)  under  a  contract  (F33615-93-C-3202),  “Smart  Materials  and  Structures  - 
Smart  Wing”  sponsored  by  DARPA. 


Figure  1-2.  Comparison  of  Pressure  Distribution  between  Hinged  and  Hingeless 

Contoured  Control  Surfaces 
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The  hinged  flap  data  depict  flow  separation  at  the  trailing  edge,  whereas  the  contoured 
flap  shows  large  suction  pressure  on  the  flap,  as  well  as  on  most  of  the  upper  surface. 
This  type  of  pressure  distribution  has  some  beneficial  effects  on  aeroelastic  stability. 
Since  the  elastic  axis  lies  downstream  of  the  section  aerodynamic  center,  the  increased 
load  near  the  leading  edge  twists  the  wing  upwards,  resulting  in  an  increased  angle  of 
attack  relative  to  that  observed  in  the  case  of  a  conventional  trailing  edge  flap.  Thus, 
the  effectiveness  of  the  trailing  edge  control  surface  increases  with  increased  dynamic 
pressure,  and  consequently,  the  roll  reversal  speed  increases.  This  is  a  significant 
contrast  to  the  case  of  traditional  control  surfaces  where  the  agility  of  the  aircraft  is 
reduced  with  increasing  dynamic  pressure. 

In  the  case  of  wind  tunnel  tests,  two  models  were  constructed.  One  model  had 
conventional  hinged  trailing  edge  flaps  and  ailerons,  while  the  other  (also  known  as  the 
Smart  Wing)  had  deformable  control  surfaces  made  of  shape  memory  alloys  (SMA). 
Typical  wind  tunnel  results  are  presented  in  Figures  1-3  through  1-5.  The  trends  are 
similar  to  those  observed  in  the  CFD  solutions.  However,  due  to  lack  of 
instrumentation,  the  pressure  loop  for  the  hinged  control  surface  was  not  observed  in 
the  wind  tunnel  data. 

Aerodynamic  coefficients  (both  lift  and  rolling  moment)  due  to  aileron  deflection  were 
measured  in  the  wind  tunnel.  The  measured  data  was  used  to  compute  the  roll  reversal 
speed  shown  in  Figure  1-6.  The  roll  reversal  speed  is  the  velocity  at  which  the  aileron 
effectiveness  is  equal  to  zero.  As  expected,  the  hingeless  control  surface  yields  higher 
roll  rate  and  higher  control  reversal  speed  compared  to  the  hinged  control  surface. 


Figure  1-3.  Pressure  Coefficient  (CP )  Comparison,  at  36%  Span 
Q  =  60  psf,  a  =  8°,  Flap  =  10°,  Aileron  =  0°(  Run  100  vs.  80) 
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Figure  1-4.  Pressure  Coefficient  (CP )  Comparison,  at  50%  Span 
Q  =  60  psf,  a  =  8°,  Flap  =  10°,  Aileron  =  0°(  Run  100  vs.  80) 


Figure  1-5.  Pressure  Coefficient  (CP )  Comparison,  at  80%  Span 
Q  =  120  psf,  a  =  8°,  Flap  =  0°,  Aileron  =  0°(  Run  109  vs.  20) 


Figure  1-6.  Aileron  effectiveness  Vs.  Nondimensional  Dynamic  Pressure 

Thus,  hingeless  control  surfaces  offer  many  improvements  over  the  traditional  control 
surfaces.  Chief  among  these  improvements  is  the  improved  aerodynamic  performance. 
Hence,  there  is  a  compelling  reason  to  investigate  various  avenues,  in  the  light  of 
smart  structures  technology,  to  develop  feasible  mechanisms  to  control  the  camber  of 
lifting  surfaces  as  is  desired  in  any  combat  mission. 

Today,  with  the  advent  of  new  materials  technology,  it  is  possible  to  design  smoothly 
deforming  lifting  surfaces  using  composite  materials.  Solid  state  actuators,  which  can 
output  large  forces  at  rates  used  in  modern  flight  control  algorithms,  are  also  being 
developed.  These  actuators  can  be  built  lightweight  in  relatively  small  sizes.  Since 
these  actuators  are  small  and  light,  a  large  number  of  these  actuators  can  be  used  on 
the  lifting  surface  so  that  any  desired  lifting  surface  deformation  shape  can  be 
commanded  for  any  flight  mission. 

To  operate  these  solid  state  actuators  with  minimum  power,  there  is  a  need  to 
determine  appropriate  actuator  locations  and  power  ratings  of  the  individual  actuators. 
This  report  discusses  a  mathematical  approach  based  on  optimal  control  theory. 
Aircraft  performance  goals,  such  as  pitch,  roll  and  yaw  rates,  are  used  as  the  target 
quantities  (constraints),  while  actuator  power  rating  is  taken  as  the  objective  function  of 
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the  design  problem.  Detailed  discussions  of  synergetic  design  methodology,  including 
the  balancing  of  the  aircraft,  are  presented  in  References  5  through  10. 


2.  ANALYTICAL  MODELING  OF  AIRCRAFT  DYNAMICS 


2.1  Overview 

The  main  objective  of  this  study  is  to  develop  an  analytical  algorithm  which  helps  to 
select  the  placement  of  actuators  and  their  power  ratings,  so  that  desired  flight  load 
conditions  and  aircraft  angular  rates  (pitch,  roll  and  yaw)  can  be  achieved  with  reduced 
actuator  weight  and  minimum  power.  A  general  approach  is  used  to  develop  a 
mathematical  algorithm  so  that  more  complex  maneuver  performance  requirements  and 
structural  design  requirements  can  be  included  at  a  later  time.  Two  mathematical 
models  are  presented,  one  based  on  the  modal  approach  and  the  other  using  reduced 
structural  degrees  of  freedom  generally  known  as  the  ASET  degrees  of  freedom.  The 
modal  approach  is  convenient  to  conduct  some  trade  studies  to  assess  system 
characteristics,  but  it  is  not  a  desirable  model  for  the  study  of  flight  dynamics  since  the 
deformation  of  the  wing  caused  by  an  actuator  cannot  accurately  be  represented  by  the 
linear  combination  of  modal  functions.  Hence,  a  brief  discussion  of  the  ASET  approach 
is  also  described. 


2.2Aircraft  Dynamics  using  Modal  Coordinate  System 

The  equations  of  motion  of  maneuvering  flexible  aircraft  can  be  written  as: 

Kr  +  Cr  +  Mr  +  QTt  Asa  +  QTt  Au  (/ K)ra  +  Fuu  +  Fb  =  0  (2-1 ) 

where: 

a  =  A  vector  of  angle  of  attack  at  aerodynamic  panels 

As  =  Steady  aerodynamic  Influence  Coefficients  (AIC)  matrix  with  respect 
to  aerodynamic  panels 

Fu  =  Nodal  forces  generated  from  the  actuator  elements 
K  =  Structural  stiffness  matrix  in  structural  degrees  of  freedom  (DOF) 

M  =  Inertia  matrix 

ra  =  A  displacement  vector  at  aero  panel  centers 

r  =  A  displacement  vector  at  structural  dof 
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T  =  A  transformation  matrix  relating  structural  dof  to  aerodynamic  dof 
u  =  Vector  of  actuator  stimuli  (input  in  volt) 

Au  =  Incremental  Unsteady  aerodynamic  influence  coefficient  matrix  for 
Fb  =  A  body  force  vector  due  to  gravitational  and  centripetal  accelerations 

^  “  -pV2  dynamic  pressure 

2 

C  =  A  damping  matrix 
k  =  col 

—  reduced  frequency,  k  >  0 

The  displacement  vector,  r,  can  be  expressed  as  a  linear  combination  of  rigid  body 
modes  and  vibration  modes.  Thus,  we  have: 


il 

EH 

n 

(2-2a) 

ra  =TTr  =  TT\in]  =  y/ar] 

(2-2b) 

where: 

yrr  =  A  matrix  of  rigid  body  modes  such  as  plunge  (Z),  sideslip  (Y),  pitch  ( 
0),  roll  (<(>),  and  yaw  (\|/).  The  control  surface  rotation  modes  may  also 
be  adjoined  to  represent  conventional  aerodynamic  effectors  if 
desired 

\ffe  =  Matrix  of  vibration  modes  (including  symmetric  and  antisymmetric 
components) 

y/a  =  Transformation  matrix  in  aero  dof 

V  ~  {tlr.Tie}  a  vector  of  generalized  coordinates, 

The  subscript  Y  denotes  rigid  body  modes,  while  ‘e’  denotes  the  elastic  vibration 
modes. 

Conventional  aircraft  use  control  surfaces  to  perform  maneuvers.  The  vibration  modal 
vectors  may  include  rigid  body  deflection  of  the  conventional  control  surfaces. 

The  angle  of  attack  (incidence)  at  the  center  of  each  aerodynamic  panel  is  given  by: 
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(2-3a) 


a  =  T—  =  —T[—  +  i/.Vr] 
V  V  ot 


-T[yrii  +  u 


dip 

ox 


V] 


or: 


where: 


a  =y/fi  +  y/xTi 

(2-3b) 

W  =  ^Ty/, 

(2-3c) 

K. 

II 

(2-3d) 

and: 

V  =  is  the  free  stream  velocity 

u  =  is  the  chord-wise  component  of  V 

%  =  {Z,Y,<p,9,ys,5p,Sq,Sr}  is  a  vector  of  generalized  rigid  body 

coordinates 

tir  =  {Z,Y,p,q,r,5q,8p,8r} ,  is  the  corresponding  velocity  vector  in 

rigid  body  motion 


The  control  surface  rotations  are  denoted  by  8.  The  rigid  body  motion  in  the  flight 
direction  (x)  has  been  omitted  since  the  computation  of  drag  is  not  accurate  in  the 
linear  aerodynamic  methods.  Should  accurate  drag  computation  tools  be  available, 
then  all  six  degrees  of  freedom  in  rigid  body  motion  must  be  included  in  the  analytical 
model. 


From  equation  (2-3b),  the  rigid  body  angle  of  incidence  and  sideslip  are  given  by: 


or,  =0-— 

'  V 

(2-4a) 

’ton  |  ^ 

1 

II 

(2-4b) 

Substituting  equations  (2-2a)  and  (2-3b)  in  (2-1),  and  pre-multiplying  by  y/T ,  one  can 
rewrite  the  equations  of  motion  in  terms  of  the  generalized  coordinates,  rj: 
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(2-5) 


Krj  +  Cfi  +  Mfj  +  QAslTj  +  QAS2ri  +  QAu(jc)rj  +  Fuu  +  Fb=0 
in  which  the  generalized  matrices  are  denoted  by  an  overbar  symbol. 

We  note  that  the  air  loads  due  to  oscillatory  motion  are  expressed  in  the  frequency 
domain.  Hence,  there  is  a  need  to  transform  the  air  loads  from  the  frequency  domain  to 
the  time  domain.  The  unsteady  aerodynamic  matrix,  Au  can  be  expressed  in  terms  of  a 
truncated  series  as: 


—  A)  ^1  +  ^  -^2  M3 

p  +  s 


(2-6) 


where: 

s  =  i  (— )  CD=iK  (2-7) 

V 

and  (3  is  the  pole  of  the  aerodynamic  lag  term. 

Using  equation  (2-6)  in  (2-5),  the  dynamics  of  the  aircraft  reduces  to: 

K'r]  +  C,fi  +  M’ii  +  QA3Xa  +Fuu  +  Fb=  0  (2-8) 

where: 

ZT=  K  +  QAsl  +  QAq 
C'=  C  +  QAs2  +  Q(^)A, 

M’=M  +  Q(-L)2A2 


(2-9) 

(2-10) 

(2-11) 


And  the  aerodynamic  lag  coordinate  is  given  by: 


X 


a 


S 

fi  +  s 


rj 


(2-12) 


2.2. 1  State  Space  Formulation 

The  second  order  differential  equation  (2-8)  can  be  rewritten  in  terms  of  the  first  order 
one  as  shown: 
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Where  state  vectors  are  given  by: 


n 


*=  n  . 


I* 


a) 


) 


and  state  matrices  are  defined  as: 


(2-14) 


(2-15) 


A  = 


0 

MrlK' 

0 


-/ 

M^C 

-I 


0 

QM  ^  A3 
V 

(y)/» 


B  = 


0 

M^FU 

0 


Fb  = 


0 

M^Fb 

0 


in  which  I  denotes  a  unit  matrix. 


(2-16) 


(2-17) 


(2-18) 


Equation  (2-13)  represents  a  general  formulation  of  aircraft  dynamics.  At  the  low  end  of 
the  frequency  spectrum  it  denotes  aircraft  flight  maneuver  problems,  while  at  the  higher 
end  of  the  frequency  spectrum  it  represents  an  aeroservoelastic  environment.  In 
subsequent  sections  the  optimal  control  theory  is  described  to  address  various  types  of 
problems  such  as: 

•  Flight  dynamics  and  flight  control 

•  Flutter  envelop  prediction  and  flutter  suppression 

•  T ransient  response  due  to  gust  encounter  and  buffeting 


2.2.2  Pade  Approximation  of  Unsteady  Aerodynamic  Matrix 
Let  the  incremental  data  be  given  by: 

4(x)  =  A(x)-A(0) 

Then  the  expansion  can  be  written  in  the  matrix  notation: 

[*][P]  =  A 

where: 

I  0 
0  I 
I  0 

R  = 

0  I 

I  0 
0  / 


(2-19) 

(2-20) 

(2-21) 


(2-22) 

(2-23) 


For  known  values  of  A(k)  at  a  number  of  reduced  frequencies,  the  AIC  matrices,  \ 
through  4  can  be  computed  from  equation  (2-20)  by  the  method  of  least  squares. 
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2.2.3  Optimal  Control  Design 


For  a  given  set  of  actuator  input  u,  equation  (2-13)  can  be  solved  for  the  response 
X  using  any  known  numerical  integration  method.  However  to  perform  flight  maneuver 
at  a  desired  flight  load  factor,  Nz,  the  actuator  input  u  will  be  computed  by  the  method 
of  optimal  control  theory  such  that  the  actuator  power  is  minimum. 

The  Hamiltonian  function  of  this  problem  can  be  stated  as: 

h=^eTQ£+U(Z)  +  AT(AX  +  BJ;  +  Fb)  (2_24^ 

where: 

A,B,Fb  =  are  defined  in  equation  (2-13) 

|  =  u  a  vector  of  control  variables 

£  =  a  vector  of  constraint  function 

Q  =  a  weighting  matrix 

A,  =  A  vector  of  Lagrangian  coefficients 

The  constraints  can  be  expressed  in  several  forms.  In  terms  of  state  variables  one  can 
write: 


€  =  (Xd-X)  =  0  (2-25a) 

Where,  in  steady  symmetric  pull-up  maneuver  at  a  load  factor,  N zr,  the  required 
constraint  is: 

%=%-<?  (2-25b) 

From  the  second  row  of  equation  (2-13)  the  nondimensional  pitch  rate  is  given  by: 

q  =  -B2t-f2b  (2-25c) 

and  the  target  pitch  rate  is: 

„  _(A^-1  )g  (2-25d) 
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Constraint  on  actuator  strain  is  given  in  the  following  expressions.  Let  ea  and  e  be 
the  allowable  strain  in  the  actuator  material  and  the  actual  strain  in  the  actuator  when 
commanded.  Then  the  difference  between  the  allowable  strain ea,  and  the  actual  strain 
in  the  actuator  must  be  positive: 

A  £  =  £a-£  =  £a-c%>0  (2-26a) 

in  which  c  is  a  constant  relating  the  actuator  strain  to  the  stimuli  (electrical  potential  in 
volts). 

The  constraints  based  on  aeroelastic  stability  criteria,  using  Lyapunov’s  direct  method, 
may  be  written  as: 

e(k)  =  ev  -  V  >  0  (2-26b) 

£  =  £y-V<  0  (2-260) 


where  V  and  V  denote  Lyapunov  energy  function  and  its  rate  of  change  in  time, 
in  which  ev  and  are  prescribed  constants. 


The  second  term  in  the  Hamiltonian  function,  eqn.(2-24),  represents  an  objective 
function  based  on  the  strain  energy  in  the  actuator  elements  which  is  given  by: 


l/(|)  =  |  (ctf  EAL(c4)  =  ^tR% 


(2-27a) 


where: 

E  =  Young’s  modulus  of  the  actuator  material 

A  =  cross  sectional  area  of  the  actuator  stack 

L  =  length  of  the  actuator 

The  weighting  coefficient  R  is  given  by: 

R  =  c2EAL  (2-27b) 

The  main  purpose  of  this  objective  function  is  to  minimize  the  power  input  so  that  the 
actuator  is  not  overly  strained  beyond  its  elastic  limit.  Likewise,  stress  or  strain 
constraints  based  on  primary  structural  members  may  also  be  included. 

The  constraint  functions  in  general  can  be  expressed  as: 
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For  symmetric  pull-up  in  steady  maneuver  we  have: 


dL  =  Q  (2-28b) 

3X 

for  all  elements  ofX , 

and  for  the  control  variables,  from  second  row  of  matrix  B ,  the  derivatives  are  given 
by: 

de-  (2-28C) 

and: 

oAe.  (2-28c) 

Additional  objective  functions  such  as  structural  weight  and  aerodynamic  figures  of 
merit,  may  also  be  included  in  the  objective  function  U . 


2.2.4  Hamiltonian  Equations  of  Motion 

Differentiating  the  Hamiltonian  function,  h  (equation  2-24),  with  respect  to  X  ,  £(=u), 

and  X ,  and  using  the  principle  of  optimal  control  theory,  we  obtain  the  following  two- 
point  boundary  value  problem: 


together  with  the  incremental  control  (design)  input: 

%  =  -*-\£?xX  +  BrX+£s  0] 
in  which: 


(2-29) 


(2-30) 


(2-31) 


The  elements  of  the  Hamiltonian  matrix,  H,  in  equation  (2-29)  are  given  by: 


(2-32) 


/in  =  A-B%  xe?c 

\2  =  -B^~XBT  (2-33) 

h21  =  —£xx  "*■  ^  (2-34) 

\i  ~  ~Ki  (2-35) 

F^-591%  (2-36) 

F2  =  —  £ xo  +  £^9^  £|0  (2-37) 

The  following  definitions  are  used: 

£v  =£iTQ£j  i  =  X,£  and  y  =  X,^  (2-38) 

£,0=^0  i  =  X,$  (2-39) 


2.2.5  Solution  to  Hamiltonian  Equations  of  Motion 

The  solution  of  the  two-point  boundary  value  problem  stated  in  equation  (2-29)  can  be 
written  as, 


f*M 

U('/)J 


(2-39a) 


where  tf  denotes  the  final  time  frame,  and  <l>  denotes  the  transition  matrix. 


The  transition  matrix  can  be  computed  as  follows. 

Let  x  and  X  be  ri9ht  and  left  eigenvectors  of  the  Hamiltonian  matrix  H ,  and  A  is  a 
diagonal  matrix  of  corresponding  eigenroots.  Then  the  transition  matrix  is  given  by: 

®(tf , t)  =  x exp(A (tf  -  t))x'  =  eH{,f~‘)  (2-39b) 

The  known  end  conditions  are: 

X(t)  =  X(t0)  (2-40) 


at  t  =  t0,  the  initial  condition,  while  at  the  terminal  condition: 
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Mtf)= o 


(2-41) 


Rewriting  equation  (2-39a)  in  matrix  notation  and  using  equations  (2-40)  and  (2-41), 
one  obtains: 

r/,(oi  <2-42> 

[_4>21  <I>22jU(Oj+l/2«.’ 

Multiplying  equation  (2-42)  by  the  inverse  of  <J>  and  rearranging  we  obtain: 

;xWi  r*»  4'12]|x(f/)i_r71(oi  (2-43) 

U(oj  K  1/2(0.' 

At  t  =  0,  we  have: 

X(t)  =  X(0) 

and  at  t  =  tf ,  we  have: 

Mtf)  =  0 


Thus,  we  get: 


(2-44) 

and  the  Lagrangian  coefficient  vector  is  given  by: 

A(f)  =  ^21'P-1{X(0}+/1}-72 

(2-45) 

or: 

A(f)  =  P(0*(0  +  J(0 

(2-46) 

in  which  the  Riccati  matrix  is  denoted  by: 

m = 'p21'p1-,1 

(2-47) 

Finally  the  control  law  is  given  by: 

m  =  -35-’[[e@  +  BrP]X(0  +  BTs(t)  +  £{0] 

(2-48) 
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2.2.6  Alternate  Riccati  Matrix  Solution 
Differentiating  equation  (2-46)  with  respect  to  time  we  obtain: 

A  =  KX  +  £X+i  (2-49) 

Next,  eliminating  X  and  X ,  using  equation  (2-29),  we  obtain  the  following  Riccati 
matrix  equation: 

K  +  Khn-h22K  +  Kh,2  K-h21  =  0  (2-50a) 

together  with: 

s  +  ( Kh12  -h22)s  +  KF ;  - F2  =  0  (2-50b) 

Since  the  initial  conditions  are  not  known  for  the  elements  of  the  Riccati  matrix  P  the 
solution  must  be  derived  from  the  terminal  conditions.  If  the  coefficient  matrix,  hi} ,  are 

constant,  the  solutions  of  P  at  various  time  intervals  may  be  computed  and  stored 
which  can  be  used  in  the  aeroelastic  analysis.  Otherwise,  the  solution  given  by 
equation  (2-47)  is  reasonably  valid  for  a  regulator  problem. 


2.2.7  Closed  Loop  System 

We  set  the  feedback  control  input: 

u=  %{t)  =  -<ft-l[[£pc  +BTK]X(t)  +  BTs(t)+£f0]  (2-51) 

and  substituting  in  equation  (2-13)  we  obtain  the  closed  loop  system: 

X  =  AX+3  (2-52) 

where: 

A  =  A  -  B<R~1BtK  -  (2-53) 

3  =  Fb  -  m~1BTs(t )  -  (2-54) 

The  eigenvalues  of  the  closed  system  A  are  stable,  assuring  the  aircraft  to  be 
dynamically  stable.  However,  appropriate  gain  matrix  coefficients,  Q ,  must  be  selected 

so  that  any  unstable  roots  that  lie  in  the  right  half  of  Nyquist’s  plane  move  to  the  left- 
half,  the  stable  region.  Lyapunov  proposed  a  robust  stability  criterion  of  a  nonlinear 
system.  A  brief  outline  is  presented  next. 


17 


2.2.8  Lyapunov  Stability  Criteria 

A.  M.  Lyapunov  (also  written  as  Liapunov)  proposed  a  robust  stability  criteria  of 
nonlinear  systems  in  1892.  This  was  practically  ignored  for  many  years.  However,  a 
French  translation  of  his  work  appeared  in  1907.  Again  it  was  forgotten  for  nearly  a 
half-century.  Due  to  a  compelling  need  to  establish  the  state  of  stability  of  nonlinear 
systems,  engineers  and  scientists  began  recognizing  the  merits  of  the  Lyapunov’s 
stability  criteria.  Consequently,  a  German  translation  appeared  in  1954  and  an  English 
translation  in  1961. 

A  Lyapunov  stability  criterion  primarily  deals  with  the  energy  and  rate  of  energy 
dissipation  of  a  system.  Hence,  it  is  valid  for  linear  as  well  as  nonlinear  systems.  There 
are  several  ways  to  formulate  Lyapunov’s  criteria.  The  stability  criteria  proposed  here  is 
well  suited  for  complex  and  multidisciplinary  systems. 

Let  the  energy  function  be  given  by: 

V(t)  =  J~  x r  PXdr  >o  (2'55) 

in  which  P  is  a  positive  definite  Riccati  matrix  solution  of  a  closed  loop  system  defined 
equation  (2-52).  The  corresponding  Riccati  matrix  equation  is  given  by: 

P  +  ATP  +  PAT  =-Q  (2-56) 

Next  we  compute  the  rate  of  change  V  by  finite  difference: 

^  =  (2-57) 

At 

A  system  is  stable  if  for  every  finite  time  t: 

V  >  0  (2-58) 

and: 

V  <0  (2-59) 

These  two  Lyapunov  stability  criteria  can  be  expressed  in  a  single  statement: 

£l=V*V<0  (2-60) 

This  concludes  the  generalized  mathematical  discussion  of  the  dynamics  of  flexible 
aircraft.  Next  we  will  present  some  applications  in  steady  maneuver  cases. 
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2.2.9  Steady  Symmetric  Maneuver 

Since  we  are  interested  in  steady  maneuver,  the  acceleration  term  in  equation  (2-13)  is 
set  to  zero.  In  addition  unsteady  aerodynamic  loads  are  not  required  in  this  study. 
Hence,  the  governing  equation  reduces  to, 


Tt  +  ^  +  Q\y\Zu  +  fb]  =  Ar)  +  Bu  +  fb  =0 


(2-61) 


where  the  state  vector  and  matrices  are  given  by: 


II 

(2-62a) 

B  =  |fe  +  0A,]"'1[f’„] 

(2-62b) 

f„=\z+QAsThfb] 

(2-62c) 

M 

(2-62d) 

77  = 


<l 


in  which  the  flight  angle  of  attack  is  given  by  equation  (2-4a).  The  inversion  indicated  in 
equations  2-62b  and  2-62c  must  be  performed  with  caution,  since  the  elements  of 
structural  stiffness  and  aerodynamic  matrices  may  be  two  to  three  order  of  magnitude 
different. 


The  expression  for  the  actuator  input  in  the  case  of  steady  maneuver,  reduces  to  a  very 
simple  form.  This  can  be  derived,  alternatively,  from  the  Newton’s  method  as  discussed 
next. 


2.2. 10  Newton ’s  Method 

The  Hamiltonian  function  will  be  replaced  by  quadratic  functions  without  using  the 
Lagrangian  coefficients  A  (see  eqn.  2-24),  i.e. 


h=±eTQxe  +  U&  +  (AX  +  BZ  +  fb)T{AX  +  BZ  +  fb) 


(2-63) 


The  resulting  equilibrium  equation  takes  the  following  form: 


where: 


SR  =  R  +  BtB  +  BlQ(2)B2  +  c2Q 
h  -  BT  fb  -cQ  +  Bl Q(2)(qT  +  fb  (2» 


(2-65) 


(2-66) 


This  approach  has  some  merit  over  the  Lagrangian  in  that  the  Hamiltonian  function  is 
strictly  a  convex  function  leading  to  an  absolute  minimum  value  of  the  Hamiltonian. 
However,  there  may  be  many  solutions.  Valid  ones  are  admitted  when  all  constraints 
are  satisfied.  In  the  present  examples  both  methods  produced  identical  results. 


2.2. 1 1  Steady  Roll  Maneuver 

Steady  roll  maneuvers  are  performed  at  1g  load  factor.  Hence,  the  external  body  forces 
vector: 


fb=  0 

and  the  equilibrium  equation  reduces  to: 


P 

Wei 


=  -[k  +  £>a]“13£  =  -B£ 


(2-67) 


(2-68a) 


or: 

Arj  +  B£=  0 

in  which  the  nondimensional  roll  rate  is  given  by: 


(2-68b) 

(2-69) 


In  this  case  antisymmetric  vibration  modes  including  one  rigid  body  roll  mode  will  be 
selected.  Then  the  generalized  coordinate  vector  is  given  by: 


The  corresponding  constraint  is  given  by: 
l=Pt-P  =  Pt~  Bi£=0 


(2-71) 


2.3  Aircraft  Dynamics  using  Reduced  Structural  Degrees  of  Freedom 

The  equations  of  motion  of  an  elastic  aircraft  can  be  written  as: 

Kr  +  Ma  +  Q'¥TAa+Fg  +  Fuu  =  0  t3'1) 

where: 

K  =  Structural  stiffness  matrix  in  structural  (ASET)  degrees  of  freedom 
M  =  Mass  matrix  consistent  with  K 
A  =  AIC  matrix  in  aerodynamic  with  respect  to  aero  panels 

=  T ransformation  matrix  relating  structural  dof  to  aerodynamic  dof 
a  =  Vector  of  accelerations  in  a  moving  coordinate  system  fixed  to  aircraft 
a  =  Angle  of  attack  at  aero  panels 
Fg  =  Vector  of  gravity  forces  and  engine  thrust 

Fu  =  Vector  of  nodal  forces  created  by  the  solid  state  actuators 

u  =  Vector  of  actuator  stimuli  (e.g.  volt) 

r  =  Displacement  vector  including  rigid  and  elastic  d.o.f. 

The  displacement  vector  can  be  expressed  as  a  linear  combination  of  rigid  body  modes 
and  elastic  degrees  of  freedom.  Thus,  we  have: 


where: 

I e  =  A  unit  matrix 

Trs  =  A  matrix  of  rigid  body  modes  at  structural  nodes 

T  =  A  transformation  matrix  defined  by  equation  (3-2a) 

re  =  {w  v  w}  a  vector  of  elastic  displacements  measured  with  respect 
to  the  statically  determinant  supports 


In  roll  and  symmetric  maneuvers  only  w component  need  be  selected  at  each  ASET 
point,  whereas  in  banked  turn  or  rolling  pullout  maneuvers  u,  v,  and  w  components 
must  be  selected. 

(3-2b) 


p={X0,Y0,Z0,j,6,y)  (3-2o) 

is  a  rigid  body  vector  consisting  of  three  displacements  and  three  rotations  of  the 
center  of  mass  of  the  aircraft.  The  corresponding  velocities  are  given  by: 

p  =  {U0,V0,W0,p,q,r}  (3-2d) 


The  trim  angle  of  attack  and  sideslip  angle  are  given  by: 


a  =  Vo. 

u0,Po  u0 


(3-3) 


The  angle  of  attack  a  is  given  by: 


r  1  ror  _rV7  .  _  1  _  .  ore 

&  =  ¥-  =  ¥-[—  +  Wr]  =  ¥\~WrP  +  u  -^]] 
V  V  ot  V  ox 


or: 

1  —  a 

<x  =  V'xr+-y'rr 

in  which  slope  transformation  matrix  from  structure  to  aero  grids: 

W*  =  [Via  0] 


(3-5a) 


(3-5b) 


and  the  rigid  body  modal  matrix  at  aero  grids: 

Wr  =  [0  yra]  (3-5c) 

The  first  term  in  equation  (3-4)  denotes  the  slope  due  to  elastic  deformation,  while  the 
second  part  represents  the  angle  of  attack  due  to  rigid  body  velocities,  with  V  as  the 
free  stream  velocity. 

The  acceleration  vector  following  Reference  5  can  be  written  as: 

a=tf+T}  (3-6) 


22 


in  which: 


^  —  { CO . . 0)  }  is  a  3n*3n  superdiagonal  matrix  with  constant 

O) 


(3-7) 


(3-8a) 


with: 


Ct)  = 


-r  q 

0  -p 

P  0 


(3-8b) 


p,  q  and  r  are  the  rigid  body  angular  velocities  of  the  aircraft  in  the  body  fixed 
coordinate  system. 

The  first  term  in  equation  (3-6)  represents  accelerations  due  to  translation  and  angular 
motions,  while  the  second  term  denote  the  centripetal  acceleration  of  a  maneuvering 
aircraft.  The  coriolis  acceleration  terms  are  omitted  being  small  in  an  inertial  frame 
fixed  to  a  flat  surface  on  the  earth. 


The  external  force  vector  Fg ,  is  given  by: 
Ft=Ma,+f 


in  which  the  gravitational  acceleration  vector  is  given  by: 

[-  SindCos(f\ 


a 


g 


Sirup 

Cos6Cos<f>  I 


a  thrust  vector 


Pitch  angle: 


6  =a  +  i 
Angle  of  attack: 


(3-9a) 


(3-9b) 


(3-9c) 


(3-9d) 
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With  these  definitions,  equation  (3-1),  after  premultiplying  by  the  transformation 
matrix,  T ,  reduces  to: 


where: 


Kr  +  M?  +  A?  +  Fg+Fuu=0 

(3-10) 

K=TTKT  +  QTTy/TA\/fx 

(3-11) 

M=TtMT 

(3-12) 

A  =  TtMTc  +  QTTy/T  Ay/ r 

(3-13) 

F  =TtF 

8  8 

(3-14) 

F  =TtF 

u  u 

(3-15) 

The  flight  dynamics  defined  by  equation  (3-10)  can  be  rewritten  in  terms  of  the  state 
space  coordinate  vector,  X ; 


X  =AX+Bu  +  fg 


(3-16) 


where: 


X  = 


A  = 

B  = 


0  -I 
-  M~lK  -M-XA 

0 

-M'lFu_ 

0 

rM-%_ 


(3-17) 

(3-18) 

(3-19) 

(3-20) 


The  solution  to  this  problem  can  be  found  by  means  of  the  optimal  control  theory 
discussed  next. 


2.3. 1  Transient  Maneuver  Analysis 
Hamiltonian  function  can  be  stated  as: 

h  =—  eT Qe  +— uT Ru  +  AT(x-Ax-  Bu  -  f )  ^ 

where  the  first  term  represents  a  function  of  target  errors,  the  second  term  denotes  the 
actuator  power,  and  the  third  term  defines  the  aircraft  dynamics.  A  is  a  vector  of 
Lagrangian  coefficients. 

The  error  vector,  e,  may  consist  of  several  constraints  such  as  the  load  factor,  roll 
angle,  angular  velocities  or  accelerations.  The  weighting  coefficient  matrix  is  given  by 

Q- 

The  constraint  function  may  be  classified  as  those  based  on  state  variables,  X  and 
the  control  variables,  u .  The  corresponding  weighting  coefficients  will  be  denoted  by 
QxandQu. 

A  typical  constraint  function  is  given  by: 

e„  =  l  -ct4,  (3-22a) 

in  which  £  denotes  either  X  or  u. 

The  angle  of  attack  is  given  by: 

a=ZL  (3-22b) 

and  the  sideslip  angle  is  given  by: 

^  (3-22o) 

P  V 


2.3.2  Hamiltonian  Equations  of  Motion 

Differentiating  the  Hamiltonian  function,  h  (equation  3-21),  with  respect  to  X  ,  u ,  and 
A ,  and  using  the  principle  of  optimal  control  theory  we  obtain  the  following  two  point 
boundary  value  problem: 


X 

A 


/iji  h j2  |X)  fFj  1 

21  ^22.1  Aj  IfJ 


(3-23) 
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together  with  the  control  input: 
u  =  -K-'[Bt?l-QuCu} 


(3-24) 


in  which: 

*  =  R  +  Qucu2  (3-25) 

is  a  diagonal  weighting  matrix. 

The  elements  of  the  Hamiltonian  matrix  in  equation  (3-23)  are  given  by: 


hn=A 

(3-26) 

\2  =  -59T1fir 

(3-27) 

^21  =  ~QxCX 

(3-28) 

•ar 

i 

ii 

-P 

(3-29) 

(3-30) 

^2  =  QxCX 

(3-31) 

The  solution  procedure  involving  two  end  conditions  at  different  time  intervals  is 
somewhat  more  difficult,  than  the  ordinary  transient  problems  with  known  initial 
conditions.  A  brief  account  of  the  solution  procedure  with  respect  to  the  transient 
maneuver  analysis  using  distributed  solid  state  actuators  will  be  presented  next  (see 
also  Ref.  6). 


2.3.3  Steady  Maneuver  Analysis 

In  the  case  of  steady  maneuver  analysis,  the  acceleration  terms,  r,  defined  in  equation 
(3-10)  will  be  set  to  zero.  In  addition,  the  contributions  from  re  and  p  can  be  omitted 

since  they  have  zero  or  negligible  effect  on  aerodynamic  forces.  Then  the  equilibrium 
equations  can  be  simplified  with  the  following  definition  for  the  state  vector: 

(3-32) 


The  corresponding  Hamiltonian  function  is: 
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h  =  -eTQe+-uTRu  +  AT(Ax  +  Bu  +  fg) 
2  2 

with 

A  =  I  a  unit  matrix 

/,=  [if +!]"?? 

The  Hamiltonian  matrix  is: 


B  =  m~lBT 

D  =  bxcx2J’  a  diagonal  matrix 

f, 

/, 

Thus,  we  obtain  the  solution: 


where: 

Cn=[I  +  Wyl 
Cv=-lI  +  DBr1 
C12  ~  BC22 
C2l  —  —DCn 

Thus,  the  actuator  input  for  required  maneuver  is  given  by: 
u=-X-1[BtA-QuCu] 


(3-33) 

(3-34) 

(3-35) 

(3-36) 

(3-37) 

(3-38) 

(3-39) 

(3-40) 

(3-41) 

(3-42) 
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2.4  Case  1: 1-g  Roll 

Given: 

1.  Roll  rate  prequire 

2.  Selected  ASET  displacement  degrees  of  freedom  w  on  the  wing  surface  and 
u  at  the  actuator  end  points. 

3.  Apply  constraints  on: 

•  Roll  rate 

•  Displacements  w  and  u 

•  Strain  in  the  actuator 

•  Strains  at  wing  root  and  mid  section  etc. 

Solution: 

1 .  Actuator  input 

2.  Net  load  distributions,  shear,  bending  and  torsion 

3.  Deformed  shape  of  the  wing 

2.5  Case  2:  Pull-up  or  Push-over  maneuver 

Given: 

1 .  Load  Factor:  Nz 

2.  Required  pitch  Rate:  q required  =^—^~ 

3.  Roll  rate:  p-  0 

4.  Yaw  rate:  r  =  0 

5.  Pitch  angle:  6  =  a  +  y 

6.  Angle  of  attack:  a  =  — 
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7. 


Climb  angle:  y  =  Sin~1(-  —) 

W 

8.  Selected  ASET  degrees  of  freedom  u  and  w  on  the  upper  wing  surface  and  u 
at  the  actuator  end  points. 

9.  Constraints  on: 

•  Pitch  rate 

•  Displacements 

•  Strain  in  the  actuator 

•  Strains  at  wing  root  and  mid  section  etc. 

Where: 

8  =  Gravitational  acceleration  inch/sec2 
V  =  Aircraft  velocity  in  inch/sec 
Solution: 

1 .  Angle  of  attack  and  actuator  input 

2.  Net  load  distributions,  shear,  bending  and  torsion 

3.  Deformed  shape  of  the  wing 

2.6  Case  3:  Banked  turn  with  Glide 

Given: 

1.  Bank  angle  ^  and  glide  angle  } 

2.  Pitch  angle:  6  =  a  +  j 

3.  Angle  of  Attack:  a  =  — 

V0 

4.  Roll  rate:  p  =  -(-^)tm$Siny 
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5.  Pitch  rate:  q  =  (■£■)  tan (pCosfiintp 

6.  Yaw  rate:  r  =  (~)  tan  <pCosyCos<f) 

7.  Selected  ASET  degrees  of  freedom  u ,  v ,  and  w  on  the  wing  surface  and  u  at 
the  actuator  end  points. 

8.  Constraints  on: 

•  Angular  rates 

•  Displacement  w  and  u 

•  Strain  in  the  actuator 

•  Strains  at  wing  root  and  mid  section  etc. 

Where: 

8  =  Gravitational  acceleration  inch/sc2 
V  =  Aircraft  velocity  in  inch/sec 
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3.  COMPUTER  PROGRAM  IMPLEMENTATION  WITH  ASTROS 
3.1  Outline  of  the  Algorithm 

The  current  implementation  of  the  Smart  Actuation  system  is  a  two  step  process, 
illustrated  in  Figure  3-1. 


Figure  3-1 .  Overview  of  ASTROS/Smart  Actuation  System 
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3.1.1  Step  One 

In  the  first  step,  a  special  version  of  ASTROS  VERSION  1 1  creates  stiffness,  mass, 
and  aerodynamic  matrices  that  are  required  in  the  Smart  Actuation  system  computer 
program.  Two  additional  bulk  data  entries,  SACNTL  and  SAELIST,  are  required  to 
define  data  that  is  specific  to  this  process.  These  bulk  data  entries  are  described  in 
subsequent  sections.  The  following  files  are  generated  by  ASTROS: 


Format 

Description 

aicmat.d 

2000(1X,E12.5) 

Matrix  of  aerodynamic  influence  coefficients. 

dyxe.d 

2000(1X,E12.5) 

ox 

Transforms  elastic  modes  to  aerodynamic 
panels. 

etype.d 

A8 

Element  types 

genkaa.d 

2000(1  X,  El  2.5) 

Generalized  stiffness  matrix. 

genmaa.d 

2000(1X,E12.5) 

Generalized  mass  matrix. 

geomsa.d 

2(18, IX), 3(1  PE12.5,1 

X) 

Steady  aerodynamics  geometry  description: 

Column  DescriDtion 

1  External  aerodynamic  box  ID 

2  Internal  aerodynamic  box  ID 

3  X  location  of  box  centroid  in  basic 
coordinates 

4  Y  location  of  box  centroid  in  basic 
coordinates 

5  Z  location  of  box  centroid  in  basic 
coordinates 

gtkg.d 

2000(1  X,  El  2.5) 

ASTROS  [GTKG]  matrix. 

parameters.d 

FORTRAN  j 

‘Parameter’ 

Problem-dependent  variable  values  which 
should  be  included  in  Kari  Appa’s  Smart 
Actuation  program. 

phi.d 

2000(1X,E12.5) 

Matrix  of  mode  shapes. 

phigt.d 

2000(1  X,  El  2.5) 

Transpose  of  ASTROS  [PHIG]  matrix. 
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psi3.d 

2000(1  X,  El  2.5) 

Rigid  body  plunge  mode  matrix  for 
aerodynamics  model.  This  file  will  only  be 
generated  for  runs  with  symmetric  boundary 
conditions. 

psi5.d 

2000(1X,E12.5) 

Rigid  body  pitch  mode  matrix  for 
aerodynamics  model.  This  file  will  only  be 
generated  for  runs  with  symmetric  boundary 
conditions. 

psir.d 

2000(1X,E12.5) 

Rigid  body  roll  mode  matrix  for  aerodynamics 
model.  This  file  will  only  be  generated  for  runs 
with  antisymmetric  boundary  conditions. 

qhhl.d 

IX, 18 

IX,  18 

2000(1X,E12.5) 

2000(1X,1PD12.5) 

Number  of  modes 

Number  of  reduced  frequencies 

Reduced  frequencies 

ASTROS  [QHHL]  matrix  list  of  generalized 
unsteady  aerodynamic  coefficients. 

This  file  will  only  be  output  if  a  flutter  analysis 
is  performed  in  the  same  boundary  condition 
as  the  SAERO  and  MODES  analyses. 

xyza.d 

18,14(1  PE12.5) 

Actuator  element  descriptions: 
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Column 

1 


Description 
Element  ID 

X  Coordinate  of  grid  1  in  the  basic 
system 

Y  Coordinate  of  grid  1  in  the  basic 
system 

Z  Coordinate  of  grid  1  in  the  basic 
system 

X  Coordinate  of  grid  2  in  the  basic 
system 

Y  Coordinate  of  grid  2  in  the  basic 
system 

Z  Coordinate  of  grid  2  in  the  basic 
system 

X  Coordinate  of  grid  3  in  the  basic 
system 

Y  Coordinate  of  grid  3  in  the  basic 
system 

Z  Coordinate  of  grid  3  in  the  basic 
system 

X  Coordinate  of  grid  4  in  the  basic 
system 

Y  Coordinate  of  grid  4  in  the  basic 
system 

Z  Coordinate  of  grid  4  in  the  basic 
system 

Young’s  Modulus 

Cross  sectional  area  or  membrane 
thickness 


3.1.2  Step  Two 

The  second  step  of  the  process  is  to  run  the  Smart  Actuation  program.  The  data  sets 
that  were  generated  in  Step  One  are  used  as  inputs  to  the  Smart  Actuation  program. 
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3.2  User’s  Guide 


3.2.1  Overview 

This  section  provides  updates  to  the  ASTROS  VERSION  11  User’s  Manual.  The 
sections  that  are  affected  are  noted  where  appropriate. 


3.2.2  MAPOL  Engineering  Modules 

This  section  provides  additional  documentation  to  Section  2.4.2. 1  of  the  ASTROS 
VERSION  11  User’s  Manual. 


MODULE 

TYPE 

DESCRIPTION 

SMARTAC 

T 

ENGINEERING 

Compute  and  assemble  data  required  for 
Smart  Actuation  calculations. 

3.2.3  Bulk  Data  Descriptions 

This  section  provides  additional  documentation  to  Section  4.7  of  the  ASTROS 
VERSION  11  User’s  Manual. 
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Input  Data  Entry  SACNTL 


Description:  Defines  control  parameters  to  the  Smart  Actuation  module. 

Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SACNTL 

SID 

ITRMAX 

RMODE 

S 

ERRTOL 

SACNTL 

101 

10 

1 

3 

0.001 

Field 

Contents 

SID  Set  identification  number  (Integer  >  0) 

ITRMAX  Maximum  number  of  design  iterations  (Integer  >  0) 


RMODES  Set  identification  number  of  a  MODELIST  bulk  data  entry  that  is  used  to 
request  the  rigid  body  modes  that  are  to  be  used  for  the  Smart 
Actuation  analysis  (Integer) 

EMODES  Set  identification  number  of  a  MODELIST  bulk  data  entry  that  is  used  to 
request  the  elastic  modes  that  are  to  be  used  for  the  Smart  Actuation 
analysis  (Integer) 

ERRTOL  Iteration  convergence  criteria  (Real,  Default  =  0.001) 
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Input  Data  Entry  SAELIST 


Description:  Defines  the  list  of  actuation  elements  to  the  Smart  Actuation 
module. 


Format  and  Example: 


1  23456789  10 


SID 

ETYPE 

BD1 

BD2 

BD3 

BD4 

BD5 

BD6 

com 

com- 

BD7 

BD8 

-etc- 

1001 

ROD 

1001 

1002 

1003 

1004 

Alternate  Form: 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SID 

ETYPE 

BD1 

THRU 

BD2 

Field  Contents 


SID  Set  identification  number.  (Integer  >  0) 

ETYPE  Character  input  identifying  the  element  type  (See  Note  1).  One  of  the 
following: 


BAR  SHEAR 

QDMEM1  TRIA3 

QUAD4  TRMEM 

ROD 
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EIDi  Element  identification  number  (Integer  >  0  or  blank) 

Remarks: 

1  •  Currently,  only  the  ROD  element  is  supported  by  the  Smart  Actuation  algorithm. 
All  other  element  types  are  included  here  for  future  development. 

2.  SID  is  here  for  future  implementations  of  the  Smart  Actuation  system  which  may 
be  more  general. 

3.  If  the  alternate  form  is  used,  EID2  must  be  greater  than  or  equal  to  EIDI . 

4.  Nonexistent  elements  may  be  referenced  and  will  result  in  no  error  message. 

5.  Any  number  of  continuations  is  allowed. 
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3.3  Programmer’s  Guide 


3.3.1  Overview 

This  section  provides  updates  to  the  ASTROS  VERSION  1 1  Programmer’s  Manual. 
The  sections  that  are  affected  are  noted  where  appropriate. 

3.3.2  Engineering  Application  Modules 

This  section  provides  additional  documentation  to  Section  5  of  the  ASTROS  VERSION 
1 1  Programmer’s  Manual. 
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Engineering  Application  Module:  SMARTACT 

Entry  Point:  SMARTACT 

Purpose: 

To  assemble  and  calculate  data  required  for  Smart  Actuation  analysis 
MAPOL  Calling  Sequence: 

CALL  SMARTACT  (  BC,  SACNTL,  SAELIST,  [PHIG(BC)],  LAMBDA, 

BGPDT(BC),  BEAMEST,  QDMM1EST,  QUAD4EST, 
RODEST,  SHEAREST,  TRIA3EST,  TRMEMEST,  [DYX], 
[GENKAA],  [GENMAA],  [AICMAT(MINDEX)], 
[AAICMAT(MINDEX)],  TRIM,  OGPWG,  CASE,  GEOMSA, 
[GTKG],  MAT1,  MODELIST,  SUPORT,  [MGG], 
[QHHLFL(BC,SUB)],  MKAEROI,  MKAER02); 


BC 

SACNTL 

SAELIST 

[PHIG(BC)] 

LAMBDA 

BGPDT(BC) 

BEAMEST 

QDMM1EST 

QUAD4EST 

RODEST 

SHEAREST 

TRIA3EST 

TRMEMEST 


Boundary  condition  number(lnteger,  Input) 

Relation  of  smart  actuation  control  parameters  (Input) 

Relation  of  smart  actuation  elements  (Input) 

Matrix  of  global  eigenvectors  from  real  eigenanalysis  (Input) 

Relation  of  real  eigenvalue  analysis  results  (Input) 

Basic  grid  point  definition  table  (Input) 

Relation  summarizing  the  CBAR  element  (Input) 

Relation  summarizing  the  CQDMEM1  element  (Input) 

Relation  summarizing  the  CQUAD4  element  (Input) 

Relation  summarizing  the  CONROD  and  CROD  element 
(Input) 

Relation  summarizing  the  CSHEAR  element  (Input) 

Relation  summarizing  the  CTRIA3  element  (Input) 

Relation  summarizing  the  CTRMEM  element  (Input) 
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[DYX] 

Elastic  modes  splined  to  the  aerodynamic  model  (Input) 

[GENKAA] 

Generalized  stiffness  (Input) 

[GENMAA] 

Generalized  mass  (Input) 

[AICMAT(MINDEX)] 

Symmetric  aerodynamic  influence  coefficients  (Input) 

[AAICMAT(MINDEX)] 

Antisymmetric  aerodynamic  influence  coefficients  (Input) 

TRIM 

Relation  containing  trim  parameters  (Input) 

OGPWG 

Relation  containing  data  from  the  grid  point  weight 
generation  computations  (Input) 

CASE 

Relation  containing  the  case  parameters  for  each  analysis 
within  each  boundary  condition  (Input) 

GEOMSA 

Relation  containing  data  on  the  geometric  location  of  the 
aerodynamic  degrees  of  freedom  (Input) 

[GTKG] 

Interpolation  matrix  relating  the  forces  at  the  aerodynamic 
degrees  of  freedom  to  the  forces  at  the  global  structural 
degree  of  freedom  (Input) 

MAT1 

Relation  containing  material  properties  (Input) 

MODELIST 

Relation  containing  lists  of  normal  modes  (Input) 

SUPORT 

Relation  containing  supported  degrees  of  freedom  (Input) 

[MGG] 

Global  mass  matrix  (Input) 

[QHHLFL(BC,SUB)] 

Matrix  list  of  generalized  unsteady  aerodynamic  coefficients 
(Input) 

MKAEROI 

Relation  containing  table  of  Mach  numbers  and  reduced 
frequencies  (Input) 

MKAER02 

Relation  containing  Mach  number  and  reduced  frequency 
pairs  (Input) 

Application  Calling  Sequence: 
None 
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Method: 


First  the  CASE  entries  associated  with  SAERO  and  SAER02  subcases  for  the  current 
boundary  condition  are  read  into  memory.  The  control  parameters  from  the  SACNTL 
entries  are  then  read.  The  number  of  normal  modes  is  obtained  by  opening  the  [PHIG] 
matrix.  The  lists  of  normal  modes  are  read  from  the  MODELIST  entries.  Then  the  list  of 
smart  actuation  elements  is  read  from  the  SAELIST  entries.  The  material  properties  are 
read  from  the  MAT1  entries.  The  Basic  Grid  Point  Definition  Table  is  read  from  the 
BGPDT  relation.  For  each  smart  actuation  element,  the  data  in  the  associated  *EST 
relation  is  read.  The  G-Set  modes  matrix  is  read,  transposed,  and  reduced  to  the 
desired  modes.  The  global  mass  matrix,  MGG,  is  read.  The  massv  vector  is  created  by 
summing  the  rows  of  MGG  and  multiplying  by  the  transpose  of  PHIG.  The  DYX  matrix 
is  read  and  reduced  to  the  desired  elastic  modes.  The  GTKG  interpolation  matrix  is 
read.  The  PHI  matrix  and  XYZA  matrices  are  assembled  for  the  smart  actuation 
elements.  PHI  describes,  for  each  mode,  the  generalized  behavior  for  each  smart 
actuation  element: 


The  XYZA  matrix  describes  the  smart  actuation  element  connectivity.  The  generalized 
mass  and  stiffness  matrices  are  read.  The  AIC  matrix  is  read.  The  free  stream  velocity 
is  read  from  the  TRIM  relation.  The  center  of  gravity  is  read  from  the  OGPWG  relation. 
The  aerodynamic  geometry  is  read  in  from  the  GEOMSA  relation.  The  SUPORT 
relation  is  read  to  obtain  information  about  the  supported  degrees  of  freedom.  For 
antisymmetric  boundary  conditions,  the  rigid  body  aerodynamic  roll  mode  matrix,  PSIR, 
is  assembled.  For  symmetric  boundary  conditions,  the  PSI3  and  PSI5  rigid  body  plunge 
and  pitch  mode  matrices  are  assembled.  If  a  flutter  analysis  has  been  performed,  the 
QHHL  matrix  list  is  read  and  reduced  to  the  desired  modes.  The  reduced  frequencies 
are  read.  The  data  required  for  the  smart  actuation  loop  is  now  assembled. 

Design  Requirements: 

None 

Error  Conditions: 


None 
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3.3.3  Database  Entity  Descriptions 

This  section  provides  additional  documentation  to  Section  9  of  the  ASTROS  VERSION 
1 1  Programmer’s  Manual. 


Entity:  SACNTL 

Entity  Type:  Relation 

Description:  Contains  control  parameters  for  the  Smart  Actuation 

System. 

Relation 

Attributes: 


NAME 

TYPE/KEY 

DESCRIPTION 

SID 

Integer  >  0 

Set  identification  number 

ITRMAX 

Integer  >  0 

Maximum  number  of  design  iterations 

RMODES 

Integer 

5et  identification  number  of  a  MODELIST  bulk  data 
entry  that  is  used  to  request  the  rigid  body  modes 
that  are  to  be  used  for  the  Smart  Actuation  analysis 

EMODES 

Integer 

Set  identification  number  of  a  MODELIST  bulk 
data  entry  that  is  used  to  request  the  elastic 
modes  that  are  to  be  used  for  the  Smart  Actuation 
analysis 

ERRTOL 

REAL 

Iteration  convergence  criteria 

Created  By: 


Module  IFP 


SAELIST 


Entity: 

Entity  Type:  Relation 

Description:  Contains  the  list  of  smart  actuation  elements. 

Relation 

Attributes: 


NAME 

TYPE/KEY 

DESCRIPTION 

SID 

Integer  >  0 

Set  identification  number 

ETYPE 

Text  (8) 

Element  Type.  One  of  the  following: 

BAR 

QDMEM1 

QUAD4 

ROD 

SHEAR 

TRIA3 

TRMEM 

EID 

Integer  >  0 

Element  Identification  Number 

Created  By: 


Module  IFP 


4.  DISCUSSION  OF  RESULTS 


4.1  Aerodynamic  Effectors 

To  verify  the  accuracy  of  the  algorithm  discussed  in  previous  sections,  a  simple  wing 
planform,  shown  in  Figure  4-1,  was  selected.  Figure  4-2  shows  the  smoothly 
deformable  trailing  edge  control  surface  mechanism  used  in  this  analysis.  Five  pairs  of 
solid  state  actuators  were  used  to  deform  the  trailing  edge  control  surface.  The 
actuators  either  pull  or  push  the  stringers.  Thus,  a  small  amount  of  bending  moment  is 
applied  to  top  and  bottom  skin  surfaces  so  that  the  control  surface  can  curl  up  or  down 
according  to  the  direction  of  the  stimuli.  In  this  model,  ten  actuators  were  used  along 
the  trailing  of  the  wing. 

An  alternate  actuation  system  using  distributed  actuators  as  the  diagonal  elements  of 
the  ribs  in  the  tip  section  of  the  wing  was  also  examined  to  determine  the  relative  merits 
of  the  two  configurations  of  aerodynamic  effectors.  Figure  4-3  shows  the  location  of 
forty  actuators.  A  number  of  performance  analyses  were  conducted.  The  results  are 
presented  next. 


Figure  4-1  A  low  aspect  ratio  wing  showing  structural  elements 
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4.2  Example  1 :  Steady  Roll  Maneuver 


Altitude 

Sea  level 

Mach  Number 

M  =  0.5 

Roll  rate 

Pt  =  3.0  radians/sec 

Number  of  actuators 

10  along  wing  trailing  edge 

Objective  function 

Minimum  total  power  required 

Constraints 

Target  value  of  pi ,  and  strain  allowable  in  the 
actuators 

The  deformed  shape  of  the  wing  in  1 .0  g  roll  maneuver  is  shown  in  Figure  4-4. 


Figure  4-4  Deformed  shape  of  the  wing  in  3.0  rad/sec  roll  maneuver 

The  integrated  air  loads,  in  terms  of  shear  force,  bending  and  torsion  moments,  are 
shown  in  Figure  4-5  through  Figure  4-7,  respectively.  The  air  load  distributions  due  to 
the  rolling  velocity,  p,  are  shown  by  thin  solid  lines,  while  the  net  load  distributions  are 
shown  by  heavy  solid  lines. 
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Figure  4-5  Shear  load  distribution  along  an  elastic  axis  located 
at  40%  chord.  M=0.5  p=3.0  rad/sec  roll  maneuver  at  Sea  level 
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Figure  4-6  Bending  moment  and  torsion  distribution  along  an  elastic  axis 
located  at  40%  chord.  M=0.5  p=3.0  rad/sec  at  Sea  level 
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Wing  Torsion  Distribution 
in  3.0  radfeec  roil  maneuver,  EA  Q  40% 


Figure  4-7  Torsion  Moment  distribution  along  an  elastic  axis  located 
at  40%  chord.  M=0.5  p=3.0  rad/sec  at  Sea  level 

The  inchworm  type  solid  state  actuators  shown  in  Figure  4-2  are  able  to  magnify  the 
displacements  by  an  order  of  magnitude.  This  is  known  as  the  mechanical  advantage 
(MA).  This  is  necessary  to  accommodate  large  relative  displacements  between  the  end 
points  or  the  attachment  points  of  the  actuators.  With  this  mechanism,  relatively  small- 
distributed  forces  will  be  able  to  deform  the  lifting  surface  so  as  to  perform  required 
flight  maneuvers.  Table  4-1  through  4-3  show  the  effect  of  the  mechanical  advantage 
on  actuator  power  requirement  and  the  net  energy. 


Figure  4-8  shows  the  energy  required  to  perform  a  3.0  radian/second  roll  versus  the 
mechanical  advantage  of  an  actuator.  For  MA  less  than  100,  a  large  amount  of  input 
power  is  required  to  stretch  the  actuator  to  match  the  relative  displacement  between 
the  end  points  of  the  actuator.  Otherwise,  desired  roll  rate  cannot  be  achieved.  For 
actuators  with  MA  greater  than  400,  the  required  energy  is  constant  at  250  ft. lbs.  This 
means  the  actuator  is  able  to  freely  accommodate  the  expansion  between  the 
attachment  points. 
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Mechanical  Advantage  =  100.0 

Net  Energy  =  1017.0  ft.lbs. 

Actuator 

Element 

Stimuli 

(Volt) 

Stress  (psi) 

Force  (lbs.) 

Energy 

(in.lbs.) 

Displaceme 
nt  (in.) 

1 

22.20 

-31,700.0 

-1,590.0 

40.0 

0.600 

2 

39,100.0 

1,960.0 

61.0 

0.280 

3 

-284,000.0 

-14,200.0 

0.495 

4 

217,000.0 

10,900.0 

1,880.0 

0.026 

5 

110.00 

-268,000.0 

-13,400.0 

-0.431 

6 

-104.00 

249,000.0 

12,500.0 

-1.494 

7 

44.19 

-133,566.6 

-6,678.3 

711.0 

-1.520 

8 

-36.21 

110,627.8 

5,531.3 

488.0 

9 

29.65 

-83,604.3 

-4,180.2 

279.1 

-2.280 

10 

-24.56 

72,489.1 

3,624.5 

209.8 

-3.740 

Table  4-1  The  effect  of  mechanical  advantage  on  the  actuator  energy  requirement,  MA 
=  100 


Mechanical  Advantage  =  500.0 

Net  Energy  =  309.0  ft.lbs. 

Actuator 

Element 

Stimuli 

(Volt) 

Stress  (psi) 

Force  (lbs.) 

Energy 

(in.lbs.) 

Displaceme 
nt  (in.) 

1 

5.200 

-10,629.8 

-531.50 

4.50 

0.598 

2 

-3.700 

788.00 

9.89 

0.280 

3 

31.170 

KfetaflIKIIl 

-6,620.90 

0.489 

4 

-21.570 

103,784.2 

5,189.20 

428.90 

1.941 

5 

22.800 

-158,524.4 

-7,926.20 

1 ,002.20 

6 

-19.900 

143,892.8 

7,194.64 

825.70 

-1.498 

7 

9.600 

-90,266.0 

-4,513.30 

324.94 

-1.514 

8 

-6.400 

73,313.2 

3,665.70 

214.34 

-2.921 

9 

6.700 

-54,929.37 

-2,746.50 

120.50 

-2.262 

10 

-4.096 

46,892.9 

2,344.64 

87.79 

-3.730 

Table  4-2  The  effect  of  mechanical  advantage  on  the  actuator  energy  requirement,  MA 
=  500 
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Mechanical  Advantage  =  1000.0 

Net  Energy  =  251.0  ft.lbs. 

Actuator 

Element 

Stimuli 

(Volt) 

Stress  (psi) 

Force  (lbs.) 

Energy 

(in.lbs.) 

Displaceme 
nt  (in.) 

1 

3.100 

-8,150.1 

2.65 

0.590 

2 

-1.330 

12,911.8 

| 

6.64 

0.270 

3 

16.030 

113,551.6 

5,677.6 

513.40 

0.480 

4 

BQgQIyH 

89,644.0 

4,482.2 

319.90 

0.011 

5 

-144,634.7 

-7,231.7 

834.30 

-0.440 

6 

hbes h 

130,595.9 

6,529.8 

680.20 

-1.500 

7 

5.290 

-84,770.3 

-4,238.5 

286.60 

-1.510 

8 

68,567.6 

MW* 

9 

-51,337.5 

105.20 

gg^M 

10 

-1.538 

43,671.5 

76.14 

-3.710 

Table  4-3  The  effect  of  mechanical  advantage  On  the  actuator  energy  requirement,  MA 
=  1000 
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Figure  4-8  Actuator  Energy  vs.  Mechanical  Advantage. 
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4.3  Example  2:  Symmetric  Pull-up  Maneuver 

Two  types  of  actuation  systems  were  exercised.  One  involving  the  deformation  of 
trailing  edge  control  surface  as  depicted  in  Figure  4-2.  Another  actuation  system  was 
represented  by  a  set  of  diagonal  members  placed  along  three  rib  sections  near  the 
wing  tip  as  depicted  in  Figure  4-3.  The  relative  merits  of  these  two  configurations  are 
discussed  next. 

4.3. 1  Case  1:  Subsonic  flow  using  the  trailing  edge  control  surface 

Symmetric  maneuver  trim  analysis  of  a  typical  fighter  aircraft  has  been  performed  using 
10  piezo-inchworm  actuators  as  shown  in  Figure  4-2.  The  flight  conditions  were: 


Gross  weight  of  the  aircraft 

10978  lbs. 

Load  Factor 

6.0  g 

Mach  Number 

0.5 

Altitude 

Sea  level 

Flight  Velocity 

6698.7  in/sec 

Target  Pitch  Rate 

0.2884  rad/sec 

Center  of  Gravity 

40%  mean  chord 

Allowable  Actuator  Strain 

0.002  inch/inch 

Table  4-4  Flight  configuration  data 


A  pull-up  maneuver  analysis  was  performed  with  the  data  shown  Table  4-4.  The  results 
in  terms  of  load  distribution,  actuator  energy  and  deformed  shape  of  the  wing  are 
presented  next.  The  deformed  shape  of  the  lifting  surface  is  shown  in  Figure  4-9.  As 
expected  the  deformed  surface  is  continuously  smooth.  The  magnitude  of  the  actuator 
forces  and  energy  required  to  perform  these  flight  maneuvers  depends  very  much  on 
the  degree  of  satisfying  the  constraint  by  selecting  the  weighting  coefficient  Q .  Large 

values  of  Q  enforce  the  constraints  very  accurately,  while  smaller  values  satisfy  within 
10%  of  the  desired  values.  Consider  for  example  the  constraint  on  the  pitch  rate 

£  —  Qj'  ^  —  0 
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For  large  values  of  the  weighting  coefficient,  Q,  the  constraint,  e,  will  be  satisfied 
exactly.  The  actuator  energy,  force  and  stimuli  for  two  values  of  Q  are  presented  in 
Table  4-5. 


Constraint 

Q 

Actuator 

Energy 

(Hp-hour) 

Actuator 

Stimuli 

(Volts) 

Actuator 

Force 

(lbs.) 

Target  pitch 
rate 

qr 

Solution  Pitch 
Rate 

<Z 

Q  =  1 

E  =  4.7e- 
08 

V=  115 

F  =  66.37 

0.00258 

0.3531 

Q 

10.0e+15 

E  =  2.797 

V 

9.9e+06 

F 

5.7e+06 

0.00258 

0.002584 

Table  4-5  Effect  of  constraint  weighting  coefficients  on  power  requirement 


In  the  first  case  when  the  pitch  rate  is  not  enforced  to  match  the  required  load  factor  the 
solution  shows  that  only  lift  equilibrium  is  satisfied  while  the  pitching  moment  remains 
unbalanced  by  a  small  amount.  For  this  case  the  actuator  input  is  115  volt  and  the 
actuator  develops  66.7  lbs.  of  force.  This  is  an  extremely  small  actuator  input  for  an 
aircraft  comparable  to  FI  6  to  perform  a  6.0g  maneuver  in  symmetric  flight. 

On  the  other  hand,  to  balance  the  pitching  moment  without  a  horizontal  tail  a  large 
value  of  the  weighting  coefficient  Q  must  be  selected  to  enforce  the  solution  pitch  rate 

q  to  match  its  target  value  qT .  This  results  in  demanding  extremely  large  actuator 

stimuli  (electrical  volts)  and  actuator  forces  as  shown  in  the  second  row  of  Table  .  The 
actuators  force  the  lifting  surface  to  deform  in  an  unusual  configuration  so  that  the 
resulting  air  load  distribution  satisfies  both  lift  and  pitch  equilibrium. 


The  corresponding  shear  force,  bending  moment  and  torsion  moment  distributions  are 
shown  in  Figure  4-10  through  Figure  4-12.  The  lift  equilibrium,  as  seen  in  Figure  4-10, 
is  achieved  irrespective  of  what  value  of  Q  is  used.  While  bending  and  torsion 
moments  are  different  for  these  two  cases. 
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Figure  4-1 1  Bending  moment  distribution  along  wing  span 
in  6.0g  pull-up  maneuver 
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Figure  4-12  Torsion  moment  distribution  along  the  wing  span 
in  6.0g  pull-up  maneuver 


4.3.2  Case  2:  Supersonic  flow  using  the  trailing  edge  control  surface 

Another  example,  using  the  same  aircraft  configuration  and  flight  conditions,  was 
conducted  in  supersonic  flow  at  Mach,  M=1.2.  The  actuator  energy  requirement  and 
actuator  force  data  are  presented  in  Table  4-6.  In  the  supersonic  case  the  unbalanced 
pitching  moment  is  very  small  compared  to  that  in  the  subsonic  flow.  Because,  the 
aerodynamic  center  in  supersonic  flow  lies  downstream  of  the  quarter  chord  and  is 
closer  to  the  center  of  mass  of  the  aircraft.  However,  even  though  the  unbalanced 
pitching  moment  is  small,  fairly  large  actuator  stimuli  and  forces  are  required  to  satisfy 
the  pitch  rate  exactly.  The  corresponding  forces  are  presented  in 

Figure  4-13  through  Figure  4-15.  Hardly  any  difference  is  seen  in  these  plots  for  the 
two  cases  with  Q=1  and  Q=10e+15. 
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Constraint 

Q 

Actuator 

Energy 

(Hp-hour) 

Actuator 

Stimuli 

(Volts) 

Actuator 

Force 

(lbs.) 

Target  pitch 
rate 

CJt 

Solution  Pitch 
Rate 

q 

Q  =  1 

E  =  5.3e-08 

V=  114.0 

F  =  65.29 

0.000448 

0.005501 

Q  = 

10.0e+15 

E  =  4.44e-04 

V  = 

9.3e+04 

F  = 

6.9e+04 

0.000448 

0.0004481 

Table  4-6  Effect  of  constraint  weighting  coefficients  on  power  requirement 


Figure  4-13  Shear  force  distribution  along  the  wing  span 
in  6.0g  pull-up  maneuver,  M=1.2 
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Figure  4-14  Bending  moment  distribution  along  wing  span 
in  6.0g  pull-up  maneuver,  M-1.2 


Figure  4-15  Torsion  moment  distribution  along  wing  span 
in  6.0g  pull-up  maneuver,  M=1.2 


4.3.3  Case  3:  Supersonic  flow  using  distributed  diagonal  actuators 

In  order  to  understand  the  nature  and  merits  of  distributed  actuators  compared  to  the 
conventional  aerodynamic  effectors,  a  set  of  actuators  was  placed  diagonally  at  three 
rib  stations  near  the  wing  tip. 
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The  aircraft  configuration  was  the  same  as  in  case  2  except  that  the  diagonal  actuators 
were  used  instead  of  the  trailing  edge  effectors.  The  actuator  stimuli,  forces  and  energy 
for  two  constraint  conditions  with  Q  =  1  and  Q  =  10.0e+15  are  presented  in  Table  4-7. 
The  corresponding  load  distributions  in  terms  of  shear,  bending  moment  and  torsion 
are  shown  in  Figure  4-16  through  Figure  4-18,  respectively.  The  deformed  shape  of  the 
wing  is  depicted  in  Figure  4-19  and  Figure  4-20.  It  is  intriguing  to  compare  the  power 
and  load  distributions  in  two  cases,  trailing  edge  effector  case  and  the  diagonal 
actuation  system.  When  Q=1  the  power  required  in  both  cases  was  practically  the 
same.  However,  the  trailing  edge  effector  case  had  smaller  unbalanced  pitching 
moment.  In  other  words  it  gave  the  air  load  distribution  which  was  closer  to  the  trimmed 
condition  both  in  lift  and  pitch.  Hence,  the  air  load  distributions  in  Figures  4-13  through 
4-15  show  very  little  difference  between  two  values  of  Q. 


Constraint 

Q 

Actuator 

Energy 

(Hp-hour) 

Actuator 

Stimuli 

(Volts) 

Actuator 

Force 

(lbs.) 

Target  pitch 
rate 

C[T 

Solution 
Pitch  Rate 

q 

Q  =  1 

4.55e-09 

63.39 

62.5 

0.000448 

0.04257 

Q=10.0e+15 

4.1449e-02 

14.67e+04 

14.47e+04 

0.000448 

0.0004481 

Table  4-7  Effect  of  constraint  weighting  coefficients  on  power  requirement 


On  the  other  hand  the  diagonally  distributed  actuator  system  shows  drastically  different 
air  load  distributions  as  depicted  in  Figure  4-16,  Figure  4-17,  and  Figure  4-18.  The  lift 
distribution  moves  towards  the  wing  root  resulting  in  reduced  wing  root  bending 
moment  by  about  50%. 


The  deformed  shape  of  the  wing  corresponding  to  6.0  g  pull-up  maneuver  in 
supersonic  speed,  M=1.2  is  depicted  in  Figure  4-19  and  Figure  4-20.  Figure  4-19 
shows  the  view  looking  from  the  trailing  edge,  while  Figure  4-20  shows  the  view  from 
the  wing  tip. 
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Figure  4-16  Shear  force  distribution  along  wing  span 
in  6.0g  pull-up  maneuver,  M=1.2 


Figure  4-17  Bending  moment  distribution  along  the  wing  span 
6.0g  pull-up  maneuver,  M=1.2 
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Figure  4-19  Deformed  shape  of  the  wing  due  to  diagonal  actuation  system  looking 

from  the  trailing  edge 
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Figure  4-20  Deformed  shape  of  the  wing  due  to  diagonal  actuation  system  looking  from 

the  wing  tip 
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5.  Conclusions  and  Recommendations 


An  analytical  simulation  algorithm  based  on  the  optimal  control  theory  has  been 
developed  to  compute  flight  maneuver  loads  using  solid  state  actuators.  A  state  space 
formulation  including  unsteady  air  loads  and  closed  loop  control  laws  has  been 
derived.  Aeroelastic  analyses  and  aircraft  flight  dynamics  can  be  performed.  A  few 
flight  maneuvers  have  been  conducted  in  1  .Og  roll  and  6.0g  pull-up  both  in  subsonic 
and  supersonic  flow  conditions.  This  study  suggests  that  the  solid  state  actuators  must 
be  able  to  travel  large  distances  between  contact  points  for  minimum  actuator  power. 
Large  displacements  and  forces  are  achievable  in  piezoelectric  devices  such  as  the 
inchworm  actuators.  To  perform  pull-up  and  pushover  maneuvers  a  canard  or 
horizontal  tail  type  secondary  load  carrying  devices  are  required  so  that  actuator  power 
requirement  will  be  a  minimum. 
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